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Abstract 

We derive the atomic radiative corrections predicted by QED using 
an alternative approach that offers the advantage of physical clarity and 
transparency. The element that gives rise to these corrections is the fluc- 
tuating zero-point radiation field (zpf) of average energy Hu)/2 per mode, 
which — in contrast with QED — is taken here as a primordial real entity 
in permanent interaction with matter and responsible for its quantization. 
After briefly recalling how quantum mechanics itself emerges as a result 
of the balance between the zpf and radiation reaction, the most impor- 
tant higher-order effects of the radiative terms on the atom are studied. 
The nonrelativistic QED formulas for the lifetimes and the Lamb shift, 
as well as the corrections to the latter due to external factors that mod- 
ify the vacuum field, are thus obtained in a self-consistent approach and 
without the need to resort to second quantization to the present order of 
approximation. 

Keywords: Radiative corrections, Atomic Lamb shift, Atomic life- 
times, Zero-point field 

1 Introduction 



The random zero-point radiation field (zpf) of mean energy hu)/2 per nor- 
mal mode, taken as a real field, has been shown in a series of recent papers 
[T]-[3] to be responsible for the basic quantum properties of matter. In par- 
ticular, the usual quantum description, as afforded, e.g. by the Schrodinger 
equation, is obtained as a result of reducing the original phase-space description 
of the entire particle-ZPF system to the configuration space of the particle in 
the time-asymptotic limit, in which an energy balance is reached between radi- 
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ation reaction and the zpf and higher-order effects of the radiative terms can 
be neglected. 

This paper is concerned with the main effects on the atom of the previously 
neglected radiative terms, calculated to lowest order in a = e 2 /He. After briefly 
reviewing some consequences of the energy-balance condition and the role of 
the ZPF in fixing the atomic stationary states, an analysis of the dynamics in 
the absence of energy balance is made, leading to formulas for the radiative 
lifetimes of excited states. Further, a calculation of the contribution of the 
radiative terms to the average energy gives the nonrelativistic formula for the 
Lamb shift. Finally, a modification of the background field through the presence 
of an external field or material objects is shown to produce in general a change 
in the radiative lifetimes and a shift of the atomic energy levels. 

The nonrelativistic, spinless, electric dipole approximation is made through- 
out the paper. The correct results (i.e. those predicted by non-relativistic 
quantum electrodynamics, QEd) are obtained in all cases. With these results 
we demonstrate that the theory of stochastic electrodynamics in its present 
form P]-[3] takes us beyond quantum mechanics, to the realm of QED. Though 
most of the results for the radiative corrections derived in the present work are 
well known, their connection with the condition of energy balance between the 
ZPF and radiation reaction is not. An interesting point is that in each case ex- 
plicit formulas for the effects discussed are obtained, along with a clear physical 
picture of their meaning 

The theory used here should be clearly distinguished from what is called 
semiclassical theory (e.g. [4]). The SED approach is not an attempt to replace 
quantum physics with a classical (or semiclassical) theory; quite the contrary, it 
is an endeavour intended to give deep physical support to quantum theory by 
answering fundamental questions as, e.g., on the physical mechanism that leads 
to the quantization and the stability of the atom, or the physical cause and 
nature of the quantum fluctuations. In this approach the ZPF is seen to play 
a crucial role for atomic stability through the energy-balance condition; there 
is no quantization in the absence of this field. The already quantized atom 
continues to interact with the zpf, which leads to the radiative corrections 
here studied. Although in the quantum regime also the field satisfies quantum 
rules, as discussed in [5], [5], these are not explicitly needed for the present 
purposes. Quantized matter under the action of the ZPF is sufficient to obtain 
the nonrelativistic radiative corrections to lowest order in a. 

2 The quantum regime 
2.1 Radiationless approximation 

For clarity in the exposition, let us recall in this section the main steps leading 
to the Schrddinger equation on the basis of the existence of the ZPF. For details 
see 0, [2]. 

The motion of the particle is governed in the nonrelativistic limit by the 
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equations (we use one-dimensional notation wherever possible, for simplicity) 



x=p/m, p = f(x) + mr'x + eE(t), (1) 

where f(x) is the external force, E(t) is the electric component of the random 
ZPF in the long-wavelength approximation and mrx is the radiation reaction 
force in the Abraham-Lorentz approximation, with r = 2e 2 /3mc 3 (w 10 -23 s for 
the electron). The density R of points in the particle's phase space is determined 
by 

m R + l {iR) + 1 {f{x) + mTXl R = ~h mR (2) 

Averaging over the realizations of the field one obtains for the (mean) density in 

E 

the phase space of the particle, R(x,p,t) = Q(x,p,t), the generalized Fokker- 
Planck equation 

& d 

— Q + LQ = e 2 —D(t)Q, (3) 

with L the Liouville operator 

f 1 9 d , . .... 

L = ——p+ —(f + mrx), (4 
TO ox op 

and D the diffusion operator, given to first order in e 2 by 

e 2 D(t)Q = e 2 J* dt' E(t)E(t>f e^-*'^. (5) 

The correlation of the electric field components is related with the spectral 
energy density of the field through 



E(t)E(t') = (4tt/3) / p(w)coau(t-l/)du). (6) 



For the ZPF (the field at temperature T = 0) the spectral energy density is given 
by 

p(uj, T) T=0 = p {w) = —^-g . (7) 

To make the transition from the phase-space equation ([3]) to a description in 
configuration space, the characteristic function Q(x,z,t) = J Q(x,p,t)e lpz dp is 
introduced, so that the marginal probability density is p{x, t) = J Q(x,p, t)dp — 
Q(x,0,t). By expanding the Fourier transform of Eq. ((3]) into a power series 
around z = and separating the coefficients of z k (k = 0, 1, 2, . . .), a hierarchy 
of coupled equations for moments of p of increasing order is obtained. The first 
two are the continuity equation and the equation for the transfer of momentum, 
which with the help of the change of variables 

z± = x ± r/z (8) 
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are shown to lead, in the limit z — > (when both z + and z_ reduce to x) 
and after some approximations, to the Schrodinger equation in terms of the 
parameter r\ (to be determined below) 

and its complex conjugate, with 

p(x,t)=iP*(x)iP(x). (10) 

It is important to stress that Eq. (|10[) is an integral part of the theory; it is 
not a subsidiary postulate. This result indicates that a regime of unitary (time- 
reversible) evolution has been attained, in which the mechanical subsystem has 
acquired its quantum properties [2J, [3J. 

2.2 The meaning of h in the Schrodinger equation 

In order for ((9|) to be fully equivalent to the Schrodinger equation, the value of 
the parameter rj appearing in it must be independent of the problem, and equal 
to H/2. Although the calculation of r\ has been presented in previous work [5J, 
[7] , we briefly reproduce it here because it serves to disclose the precise point of 
entry of Planck's constant into the Schrodinger equation^] 

The value of r\ will be determined by resorting to the energy-balance con- 
dition, which equates the average power lost by the particle through Larmor 
radiation, to the average power extracted by the particle from the random field. 
To establish the energy-balance condition we take the generalized Fokker-Planck 
equation ([3J, namely 

+ -JTVQ + tt(/ + mrx)Q = e 2 —D(t)Q, (11) 
at m ox op op 

(p = mi) multiply it by p 2 and integrate over the entire particle phase space. 
Assuming all surface terms to vanish at infinity, we obtain 

1 d i o\ 1 d f n 1/ ... n ~ 

— — (p z ) = — — / p z Qdxdp= — (fp + mrpx - e>D 
2m dt x ' 2m dt J m \ 

where (g) = J g(x,p)Qdxdp. Since d(V)/dt = — (fp) /m, the total average 
energy gain or loss per unit time is given by 

l^ = l{L p2+v ) =mT{± * ) -i( pt) )< (12) 

where H is the mechanical Hamiltonian function.. The first term on the right- 
hand side represents the average power dissipated by the particle through Lar- 
mor radiation; the second term represents the average power extracted by the 



1 In addition, the present calculation serves to correct a factor 1/2 mistakenly introduced 
in previously published versions of Eq. \12\ and the following. 
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particle from the ZPF and absorbed by the momentum fluctuations. For energy 
balance to hold in the mean, these terms must compensate each other, i.e., 
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mrixx) = — (pD) . (13) 
m \ / 

For a calculation of these two average values to lowest order in t ~ e 2 we use 
the solutions of Eq. © containing the parameter 77. Since the zpf represents 
the background field in its ground state, the particle must also be in its ground 
state (denoted with the subindex 0). For the left-hand side of Eq. (IT3l) this 
leads to 

mT(xx} = -mr^w„ t \x k\ 2 , (14) 

k 

where wofe = (£0 — £fc) / 2?7, x$ k = J ipo x1 Pkdx, £fc are the energy eigenvalues 
and il) k the corresponding eigenf unctions. For the calculation of the right-hand 
side of Eq. (fT3|) . which is somewhat more elaborate, we introduce Eq. (0 for 
the spectral energy density of the field — which is proportional to h — into Eq. 
((5]), whence 

— IpD) = — [duju: 3 [ dt'coSLj(t-t')I(t-t') (15) 
m \ / 7r J J 

with 

I(t-t') = Jdx J dppe-^-^-^-Qit') = JdxJ dpp-^Q{x',p',t'), (16) 

where x' ,p' are the position and momentum variables, respectively, which evolve 
deterministically (under the action of L) towards their final values x — x(t), 

p = p{t). Upon integration by parts, and writing J dxdp = J dx' dp' to zero 

order in e 2 , we get 

/ (*- i ') = <^) =^{[x',p]) = -Y,Ukv\xok\ 2 cosu Jko {t-t'). (17) 
\ op ' / 2iri Tj ^ 

Now we insert (fl7|) into (|T5|) and integrate over time starting at —00 (with 
y = t — t'), to take into account that energy balance is established after particle 
and field have interacted for a sufficiently long time, i.e. in the so-called time- 
asymptotic limit 

2 00 

— (pi)) = ^ w fe0 |x fc| 2 / du a; 3 / dy coswycos uikoy 

m \ /a ttt] J J 

K 

00 

= — ^7"5Z a;fe0 \ x ^\ 2 J du; to 3 [5 {uj + u k0 ) + S{uj - cj fe o)] (18) 
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For Wfco > (as is the case for the ground state) the first integral is nil, whence 



sH-t? 4 ^- (19) 

On comparing with Eq. (|14p we obtain 

v = n/2, (20) 

and Eq. Q becomes the Schrodinger equation, 

"Te—is?* +v *- (21) 

Note that the ZPF has played a crucial role in leading to this result. Firstly, it is 
the source of the Planck constant in this equation, through the spectral energy 
density given by Eq. ([?])• Further, a field with energy spectrum proportional to 
w 3 (responsible for the oj^ k factor in Eq. (fT!?)) ) is the single one that guarantees 
detailed balance, by ensuring that Eqs. (TH)) and (ITO1) have exactly the same 
structure. This means that energy balance holds not only globally but term by 
term, or for each frequency. This differs essentially from the result obtained 
for a classical multiply periodic system in equilibrium with a radiation field, in 
which case balance is attained only if the field has a Rayleigh- Jeans spectrum, 
proportional to w 2 [8]. 

2.3 Detailed balance for an excited harmonic oscillator 

Let us now consider an atom in an excited state n, with the background field 
still being in its ground state (the ZPf). Then instead of ([M)) we have 

mr{xx) n = -mT^uJ t \x nk \ 2 , (22) 

k 

and instead of (flU)) (with 77 = fi/2) we have 
2 

'pf)\ = -mT^u 4 nk |x„ fc | 2 sign(a; fe „). (23) 

" k 

Whilst in Eq. (|2"2"T) all terms have the same sign, Eq. (|2"3"|) contains now a 
mixture of positive and negative terms. As a result there is a net average loss 
of energy 

dt " ~d~t~ = ~ mT X! w ™fc l x " fc | 2 (1 ~ sign(wfe„)) = -2mr^w^ \x nk \ 2 . 

k k<n 

(24) 

Hence, as was to be expected, there cannot be detailed balance between the 
zpf and the atom in an excited state; in the sole presence of the ZPF, only the 
ground state (n = 0) is in equilibrium. 
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Let us now assume that also the background field is in an excited state, and 
inquire whether in this case there can be equilibrium between particle and field. 
We write the spectral density of the excited field as 

pH = PoMtM, 7H > i 
where the additional contribution 

Pa(u) = P - Po = Po(l - 1) = Pola (25) 

can represent an excitation of the background field or an external field. Now 
observe that the generalized form of Eq. (1231) for the case 7(0;) > 1, 

— (pD) = -mr y^uJn k |x„ fc | 2 7(|a;„fc|)sign(a;fc„), (26) 



m 



k 



contains again a mixture of terms with different signs depending on the sign of 
LUkm so that 

= -tot^u;^ \x nk \ 2 [1 - 7(|w„ fc |)sign(w fen )]. (27) 

k 

When the values of \ui n k\ differ for different k, the positive and the negative 
terms in this equation cannot compensate each other in general, and detailed 
balance can therefore not be satisfied. However, if all values of \ui n k \ in Eq. (|27P 
are equal, it may be possible that detailed balance exists with the particle in an 
excited state n. This is precisely the case of the harmonic oscillator: all \ui n k\ 
appearing in (I27[) are equal and moreover coincide with the oscillator frequency 
ujq. With |x„ n +i| 2 = a(n + 1), |x„„_i| 2 = an, and \x n k\ 2 = for k 7^ n ± 1 , 
where a = fi/2mwo, Eqs. (f2"2"j) and (|2"6"|) give 

1 e 2 / - \ 1 

mr{xx) =--HTul(2n + l), — (pD) = --firwg7(wo), 

whence the energy gain or loss is given by 

^L = -ifi 71u g[(2n + l)-7M]. 

Therefore, detailed balance exists between a harmonic oscillator in its excited 
state n and an (excited) background field with 7(0;) = 2n + 1. According to Eq. 
(|23|) . this field has a spectral energy density 

pju) = Po (oj)(2n + 1) = ^3 (2n + 1), (28) 

corresponding to an energy per normal mode ifiw(2n +1), equal to the energy 
of the mechanical oscillators with which it is in equilibrium. This is simply the 
condition for balance between field and matter oscillators of the same frequency 
— and a result that links with the Planck distribution in the case of thermal 
equilibrium (see [6] and section 13721 below) . 
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3 Spontaneous and induced transitions 



3.1 Radiative lifetimes 



Equation (|27|) determines the average rate of energy loss or gain by the mechan- 
ical system in an excited state n due to (upward or downward) transitions to 
states k with k > n and k < n, respectively. In order to analyze this equation 
in detail it is convenient to write 7(0;) = 1 + 7„(w), as follows from Eq. ([25]). 
and separate the positive from the negative terms, 

= -mr^u^fc \x nk \ 2 [1 - (1 + 7 a (|w„ fc |))sign(w fen )] 

k 

= mr^u^fc \x nk \ 2 [(7aW„>o - (2 + 7aW„<o] ■ (29) 

k 

The first term within brackets in (|29p represents the upward transitions (ab- 
sorptions) and the second one, the downward transitions (emissions). It is clear 
that upward transitions can take place only when there is an additional field 
7a( w fen) from which the atom may absorb the necessary energy; in other words, 
the atom does not ('spontaneously') absorb energy from the ZPF. This is an 
important point that explains, for example, why optical detectors, including 
photographic plates, are not activated by the vacuum. Emissions, on the other 
hand, can be either 'spontaneous' (in presence of just the ZPf) or else stimulated 
by the additional field j a (tUkn), according to the second term in ([29|l . Notice in 
particular that the ground state (n = 0, uj kn > 0) is absolutely stable against 
spontaneous transitions. 

Let us now relate the coefficients appearing in the various terms of Eq. (|2^|) . 
with the respective Einstein A and B coefficients. By writing the rate of energy 
change in terms of the latter 

= y^h\u} nk \ [p a (\uj nk \)B kn ] - ^h\u nk \ [A nk + p a (\u nk \)B nk ] (30) 

k>n k<n 

and comparing term by term with (|29[) , we obtain for the 'spontaneous' emission 
coefficient 

A 2Tm I I 3 I I 2 4e2 Kfcl 3 I |2 /o-n 
Ank = — Y~ \UJnk\ \Xnk\ = ^3 \X n k\ (ol) 

and for the stimulated transition coefficients 

_ _ TUT \uj nk \ 4 |x»fc| 2 7q(l^nfc|) _ 47T 2 e 2 2 

£>nk — -Dfcn — t I in l\ — ofc? \ x nk\ 7 V" 3z V 

in full agreement with the respective QED formulas [9j.|10j. 

3.2 Spontaneous decay and the zero-point field 

One can frequently find in the literature that all the spontaneous decay is at- 
tributed to cither the vacuum fluctuations or radiation reaction, more often to 
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the latter (see e.g. ([10] . [H]-[I1]). Let us look at this issue from the perspective 
of the present theory. 

From Eqs. (I5T|) and the ratio of the A to B coefficients is 

,U ftk^ (33) 

Bnk n^c 6 

Incidentally, this relation and the equality of the coefficients B nk = B kn were 
predicted by Einstein on the basis of his statistical considerations ([15]; see 
below) . 

Notice in particular the factor 2 in equation Given the definition of the 
coefficients, one could expect the ratio in this equation to correspond exactly to 
the spectral density of the zpf, which would mean a factor of 1. However, as 
follows from Eq. (f29|) . one should actually interpret the factor 2 as 2 = (1 + 1). 
One of these two equal contributions to spontaneous decay is due to the effect 
of the fluctuations impressed on the particle by the ZPF; the other one is due 
to Larmor radiation. Their equality (with opposite signs) leads to the exact 
balance of these two contributions when the particle is in its ground state, thus 
guaranteeing its stability (cf. Eq. (IM1) for n = 0). 

A brief digression is in place regarding the point at which Einstein introduced 
quantization in his 1917 paper |15j . so as to arrive at the Planck distribution. 
It is frequently argued that he did so through the assumption of discrete atomic 
levels. However, some time after Einstein's original work, Einstein and Ehren- 
fest [IB] showed that this was not the case, by redoing the calculations with a 
continuous distribution of atomic levels. In line with the results presented here 
and in previous work pQ, [3] quantization enters through the introduction of a 
source that includes the ZPF, able to generate 'spontaneous' transitions. This 
can be easily verified by omitting in the calculation any of the three terms that 
lead to matter-field equilibrium: stimulated absorptions and emissions, or spon- 
taneous emissions. The absence of the latter leads to absurd results, as happens 
also with the omission of stimulated absorptions. The omission of the term 
related to stimulated emissions leads to the expression for the blackbody law 
proposed by Wien, which correctly approximates Planck's law at low tempera- 
tures, so it already contains some quantum principle due to the presence of the 
term associated with spontaneous emissions. All this can be easily seen in the 
present context by focusing on just two states n and k, with £ n — £ k = fu^nk > 
and respective populations 7V ra , N k . When the system is in thermal equilibrium 
at temperature T, the relation [ks is Boltzmann's constant) 

N k /N n = exp(£„ - £ k )/k B T 

holds (disregarding inconsequential degeneracies). Since according to Eq. (|29p 
the number of emissions is proportional to N n "f a (uj nk ) and the number of absorp- 
tions is proportional to N k [2 + 7 (w n fe)], from the (detailed) balance condition 
Nnla — A^fc(2 + 7 a ) one obtains indeed Planck's law (for the thermal field) as 
is well known since 1917 [15] . 

2 

lai^nk) = 77 p \ lu rr 7- ( 34 ) 

exp(£„ - £ k )jk B T - 1 
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4 Radiative corrections to the energy 



The determination of the Lamb shift has been one of the most frequently studied 
problems in SED and has produced some successful results in the past, though 
basically restricted to the linear-force problem. Early related works are |17j - 
[20] ; additional references can be seen in [5T] and [35]. The theory of SED as 
developed recently and used in this paper, has the advantage of being applicable 
to nonlinear forces in general and to the atomic problem in particular. Since this 
theory includes the radiative terms from the outset, we can use it also to derive 
general formulas for the radiative energy corrections. Here we present a full 
derivation of the nonrelativistic atomic Lamb shift and associated effects within 
the present framework. The results obtained are in line with the predictions 
deriving from QED; however, the procedure followed for their derivation offers a 
clear picture of their physical meaning and allows a comparison with alternative 
interpretations of the Lamb shift found in the literature. 

4.1 The Lamb shift 

To calculate the energy corrections we go back to the generalized Fokker-Planck 
equation (jlip and multiply it now by xp before integrating over the entire phase 
space; the result is (assuming again all surface terms to vanish at infinity) [23] 

— {xp) = — (p 2 ) + (xf) + mr (x x) - e 2 (xD) . (35) 

In the radiationless approximation, corresponding to quantum mechanics, the 
last two terms are neglected and Eq. ([35]) reduces to 

J t («P> = ^ (P 2 ) + <*/> • (36) 

Further, in a stationary state (denoted again by n) the time derivative of (xp) 
is zero and Eq. (|36[) reduces to the virial theorem, 

1 </> n = 2 <T)„ = -(*/)„, (37) 

where (T) is the average kinetic energy. Hence Eq. (|35l) can be interpreted 
as a time-dependent version of the virial theorem, with radiative corrections 
included. Observe that here the average is taken not over time, but over the full 
particle phase space, which is equivalent to an ensemble average. This is but 
an example of application of the ergodic properties acquired by the quantum 
states as discussed in detail in [5], [BJ. 

In the stationary state, the two previously neglected terms can therefore be 
taken as radiative corrections to the (kinetic) energy, 

5£ n = 6 (T) n =~^{x x) n + ^ (xt)) n . (38) 
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This is a general expression for the Lamb shift. Notice that the general laws 
derived here — such as (fT2j) or ([38]) — are foreign to quantum theory, where the 
notion of diffusion (and the related diffusion operator) does not appear at all. 
The right 'hand side of Eq. (|38p will again be calculated to lowest order in 

e 2 , which means calculating the two average values, (x x) n and ( X -D^ ; to zero 

order in a — e 2 /he. For the first one we get 

rriT t , , r d , , 

— <**>.= 2 <*/>»= 2 5 CT>» = 0, 

which means that the Larmor radiation term does not contribute to the energy 
shift in the mean, in a stationary state. The correction to the energy comes 
exclusively from the fluctuations due to the action of the background field on 
the particle, represented by the second term in Eq. (f38|). To calculate it we use 
again Eq. ([5]), multiply it this time by x, 

2h „a [ j. . f ^,,,3_^../ t _^\_--£(t-t'). fl 



e 2 xD(t)Q = —^e 2 / duo / dt'cu 3 cos u(t - t') x e~ L(t ~* > — Q(t'), 
37rc cl J J op 

and integrate over phase space, thus obtaining 

e 2 / » \ he 2 f f 

Y\ xD ) =3 — 3 / dujuj3 / dt'cosu}(t-t')J n {t-t'), (39) 

with 

J n (t-t') = J dx J dpxe-^-^-^-Qit') 

Upon an integration by parts (again with J dxdp = J dx' dp 1 to zero order 
a) we have J n (t —t') = — (dx/dp') n , where 

/ dx \ 1 2 a 

Jn{t - if) = = -ft (tM'])n = ~J^\x nk \ Sm-OJknit - t'). 



Inserting this result into p9p we obtain 

°° t 

^2 \x n k\ 2 J du uj 3 J dt' cosw(t-t')$mcJkn(t-t'). (40) 



iut>\ . &2 



2 \ In 3ttc 3 

k 



Extending the initial time integral to — oo (as corresponds to the time- asymptotic 
limit) we have (with y = t — t') 

/>oo 

/ dy cos uy sin oj kn y = (41) 
Jo 

1 f°° w fc 
' dy [sin(o; fe „ + w)y + sin(o; fc „ - u)y] = — — — t , 



2 



J kn 
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which introduced in Eq. (|4T)]) gives for the radiative correction to the (mean ki- 
netic) energy (we write the result in three dimensions, for comparison purposes) 



,3 



5£ n = -{ x -D) = -— -^|x nfc | 2 Wfe „ / duj -. (42) 



This result coincides with the formula derived by Power [53] for the Lamb 
shift on the basis of Feynman's argument [25]. We recall that according to 
Feynman, the presence of the atom creates a weak perturbation on the nearby 
field, thereby acting as a refracting medium. The effect of this perturbation is 
to change the frequencies of the background field from u to u)/n(uS), n{uj) being 
the refractive index. The shift of the ZPF energy due to the presence of the atom 
in state n is then ([TO], [24"]) 

A£„ = V ; kn ; — — fujl-n — — "V^ fan (Wfcn) _ ll — fiwtn, 

^2!i„(u h ) ^2 "V knj J 2 fcn, 

where a summation over the polarizations is included, and the refractive index 
is given in this approximation by 

+ (43) 

where d mn = ex mn is the transition dipole moment. After an integration over 
the solid angle and summation over the polarizations, Power obtains in the 
continuum limit for Uk the formula 



3 



2 2 : 



which is equal to our result, Eq. (|4"2"]l . 

The Lamb shift proper (the observable Lamb shift) is obtained by subtracting 
from the total energy shift given by Eq. (|42|). the free-particle contribution S£[ p . 
This latter is represented by (|4"2")l in the limit of continuous electron energies 
(when u>k n can be ignored compared with u in the denominator), 



SSfp = ^3 ^ |x„ fc | 2 w fcn / du) lo = y w. (44) 

fc 

To write the last equality we used the sum rule £fc |x n fc| Wfc„ = 3H/2m. Since 
according to Eqs. ^ and ([7]) 



A 2 B = 2* / cLo lo, (45) 



7TC 







12 



with A the electromagnetic potential associated with the zpf, (|4"4")l can be 
rewritten as 

S£ ip = -^A^, (46) 

which identifies this 'free-particle contribution' to the Lamb shift with the con- 
tribution from the (free) background held. Its value is independent of the state 
of the particle; it is a testimony of the ubiquitous presence of the zpf. Inserting 
the usual cutoff frequency lu c — mc 2 /H in the integral in Eq. (|4"4")l gives the 
finite result 6£f p = (a/2ir)mc 2 . 

By subtracting (l4"4")l from (l4"2l we obtain for the Lamb shift proper 

2 °° 

S£ hn = 8£ n - 5£ {p = -- — - V |x„ fe | 2 ul n / du — r . (47) 

k 

Inserting once more the cutoff frequency u> c — mc 2 /h in the integral givetH 

(48) 



k 



Huii 



which is Bethe's well-known expression [26] . 

An important difference between the procedures used in the present paper 
and in QED to arrive at the Lamb shift formula concerns the mass renormal- 
ization. We recall that in the QED case, second-order perturbation theory is 
used, with the interaction Hamiltonian given by H lnt — — (e/mc)A • p. But the 
energy derived from this term, namely jlOj 



oo 

2e 2 



o 3 / ^ nK kn I • 

k 

still contains the (linearly divergent) free-particle contribution 
2e 2 



3ttc 



5>„ fc i 2 < (^Eip^ij / 

fc \ fc / 



du 



that must be subtracted to obtain the Lamb shift proper. Because this result 
is proportional to the mean kinetic energy, the ensuing correction is taken to 
represent a mass renormalization, 



4e 2 
3irc 3 



du, (49) 



2 Note that to get a correct (and finite) result, it is essential to leave in the denominator 
of this formula the term T 2 ujf nn ~ t 2 w 4 under resonance due to the presence of radiation 
reaction (I17|, |27p. This is a natural term in both QED and SED. 
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which with the usual cutoff oj c = mc 2 /h becomes 8m = (Aa/iii)m. 

By contrast, in the derivation presented here to obtain <5£l„, Eq. (|47p . there 
was no mass renormalization. The result (|49[) is just the classical contribution 
to the mass predicted by the Abraham-Lorentz equation ([22], Eq. 3.114) (or 
Maxwell's equations). In the equations of motion ([T]) this contribution has been 
already subtracted, so there is no more need to renormalize the mass. However, 
and as is well known, the formula (|47[) (common to both SED and renormalized 
QED) still has a logarithmic divergence that can be remedied by introducing the 
cutoff frequency u) c , as was done by Bethe, thus obtaining a very satisfactory 
result for the Lamb shift. 



4.2 Alternative interpretations of the Lamb shift 

The interpretation of the Lamb shift as a change of the atomic energy levels 
due to the interaction with the surrounding zpf is fully in line with the present 
theory. It constitutes one additional manifestation of the influence of the particle 
on the field, which is then fed back on the particle. An alternative way of looking 
at this reciprocal influence is by considering the general relation between the 
atomic polarizability a and the refractive index of the medium affected by it, 
which for n(ui) ~ 1 can be written as follows, 

n(to) = 1 + 2na(uj). 

Comparing this expression with Eq. (|43l) we obtain 

6h ^ u> mn - LU 



which is the Kramers-Heisenberg formula [TT]. This indicates that the Lamb 
shift can also be viewed as a (second-order in eE) Stark shift associated with 
the dipole moment d(w) = a(w)E induced by the electric component of the zpf 
on the atom. 

Equation (|47|) can be further recast in an (approximate) form that is usual 
to find in textbooks and Lamb-shift calculations, by assuming that the inte- 
gral depends so weakly on the index k that such dependence can be ignored. 
Expressing (I47|) in terms of energy levels, with huj n k = £ n — £k, one gets then 

2 °° 

S£hn = -7r^^2\xnk\ 2 u 3 kn d£ 

I., J 



k 







(£ k - £ n Y - £ 2 



= — n 2 2 / \Pnk ■ Pkn\ \£k ~ £n) ■ 

which after a series of transformations becomes [TU] 
2al n y ah 2 I n 



37TC 2 TO 2 37TC 2 m 2 



V 2 V 
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This result suggests to interpret the Lamb shift as due to the variations of the 
potential energy originating in the fluctuations in x-space. For the Coulomb 
potential, X7 2 V = 4TrZe 2 S 3 (x), so only the wave function at the origin (s states) 
contributes to the Lamb shift in this approximation. This makes this formula 
particularly practical for numerical calculations. 



4.3 External effects on the radiative corrections 

From the results obtained above it is clear that certain basic properties of the 
vacuum — such as the intensity of its fluctuations or its spectral distribution 

— are directly reflected in the radiative corrections. This means that a change 
in such properties can in principle lead to an observable modification of these 
corrections. The background field can be altered, for instance, by raising the 
temperature of the system, by adding external radiation, or by introducing 
objects that affect the distribution of the normal modes of the field. 

Such external or 'environmental' effects have been studied for over six decades, 
normally within the framework of quantum theory, although some calculations 
have been made also within SED, again for the linear- force (or single- frequency) 
problem only, leading to comparable results (see [2E],[2H]). The results of the 
previous sections, by contrast, can be applied to the general case, without re- 
stricting the calculations to the linear-force problem. In the following we present 
an illustrative selection of results derived for both lifetimes and energy levels. 
The task is facilitated by the use of the present theory because the influence of 
the background radiation field (and its modifications) is clearly pictured from 
the beginning. 

In section 13.11 we have already come across one observable effect of a change 
in the background field: according to Eq. ([29)) the rates of stimulated atomic 
transitions are directly proportional to the spectral distribution of the external 
(or additional) field, be it a thermal field or otherwise. In the case of a thermal 
field in particular, with 7 (|w n fc|) given by Eq. (IMl) . the (induced) transition 
rate from state n to state k becomes (using Eqs. (f30f and (|32l) ) 

dN nk 4e 2 \u nk \ 3 \x nk \ 2 1 

= Po(\Unk\h a (\u>nk\)B nk = — 3 e K\^ k \/k B T _ i ( 51 ) 

This result shows that no eigenstate is stable at T > — as is well known 

- because the thermal field induces both upward and downward transitions. 
For downward transitions (uj n k > 0) we can rewrite Eq. ([5"Tj) for comparison 
purposes in terms of A nk as given by Eq. (|31|) . obtaining 

dN nk _ A nk 
dt ~ e h \^ k \/k B T _i ^ > 

At room temperature (ksT ~ .025 eV) the effect of the thermal field on the 
decay rate is barely noticeable, since for typical atomic frequencies the inverse 
of the denominator, (exp h \uj n k\ /ksT — 1) , ranges between exp(— 40) and 
exp(— 400). The decay of excited states is therefore mostly spontaneous in this 
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case. For the thermal field to have a noticeable effect on the decay rate, the 
temperature would have to be of the order of 10 4 K, at which other effects on 
the atom (assuming it still exists at this high temperature) cannot be ignored. 

When the geometry or the spectral distribution of the field is modified by 
the presence of conducting objects (such as metallic plates or the walls of a 
cavity), the transition rates are affected accordingly. Let us assume, for sim- 
plicity, that the modified field is still isotropic, with the density of modes of 
a given frequency \to n k\ reduced by a factor 7 a (|w n fc|) < 1: then according to 
the results of section [331 the corresponding spontaneous and induced transition 
rates are reduced by precisely this factor, since both A and pB are proportional 
to the density of modesH By enclosing the atoms in a high-quality cavity that 
excludes the appropriate modes one can therefore virtually inhibit the corre- 
sponding transition. For the more general (anisotropic) case the calculations 
are somewhat more complicated, without however leading to a substantial dif- 
ference from a physical point of view. Such cavity effects have been the subject 
of a large number of fine experimental tests since the early works of Kleppner 
and others within QED ([30]- [33]; for more recent work see, e.g., [54"] . [35] ) . 

For illustration purposes, let us also briefly indicate how Eqs. (|44t and (|47| 
can be used to calculate the changes in the atomic energy shifts produced by 
the addition of an (external or thermal) field. As in Eq. ([23]). we denote by 
Pa = Pala tne spectral (energy) density of the additional field. The formulas 
for the variations of the (first-order) radiative corrections are readily obtained 
by determining the shifts produced by the total field (p + p a ) and subtracting 
the original shifts produced by the ZPF alone. The results are 

oo 

A (5%) = -^ [duw, (53) 



o 



A(^ Ln ) = -^£|x„ fc | 2 u,£„ / <kj 7 U (54) 

ft Q 

for a homogeneous, isotropic field. If the additional field represents blackbody 
radiation at temperature T, 7 a is given by Eq. (|34]> . i.e. 7 a (T) = 2/(expy — 1) 
with y = (fiaj/fcsT), and we obtain from Eq. (|53|) 



oo 

A T (5£ ip ) = -^{k B Tf [ dy V —-, (55) 

nmc z J exp y — 1 

o 



whence the free-particle energy increases by the amount 



A T (6£ ip ) = —(k B T) 2 . (56) 



3 Interestingly, however, by virtue of this proportionality, the ratio of spontaneous to in- 
duced transition rates is not altered by a modification of the density of modes. 
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The formula for the change in the Lamb shift is given from Eq. f[M|) by 



oo 



A(^l„) 



4e 2 





) 



37TC 3 



duj 



These results coincide with those obtained through considerably more cumber- 
some procedures within QED |36j , |37j , and the corresponding thermal shifts have 
also been experimentally confirmed ([35]; see also From the point of view 

of SED (or QEd) their interpretation is clear: they represent additional contri- 
butions to the kinetic energy impressed on the particle by the thermal field, 
according to the discussion at the beginning of section 14.11 

5 Concluding remarks 

All results contained in this paper for a nonrelativistic spinless particle point 
to the zero-point radiation field as the source not only of quantum behavior 
itself, but also of the radiative effects on quantum systems. The stochastic 
problem posed by the action of the ZPF on the particle led to a generalized 
Fokker-Planck equation for the particle, and all results presented here have 
been directly derived from this equation. Although as a result of the matter-ZPF 
interaction both matter and field end up quantized, for the present calculations 
of the radiative corrections (to lowest order in a) it sufficed to consider the 
ZPF as a classical function. Indeed, the theory furnishes an alternative way to 
derive self-consistently results usually considered to be the exclusive province of 
QED. Important advantages of the present procedure are physical transparency 
and simplicity; there was no need to resort to heuristic arguments along the 
derivations. These advantages are particularly apparent in the calculation of 
environmental effects on the atomic lifetimes and energy levels. 

It is important to stress that the present theory implies quantization of 
both matter and radiation field [I], [2J, [5J. This is the ultimate reason that 
guarantees its equivalence with QED. The difference between these two theories 
lies not in their final results, but in the whole conceptual picture and the gained 
clarification of the physics. The present approach gives well-defined answers 
to deep questions, such as the origin and ultimate meaning of the Schrodinger 
equation and other puzzles of quantum theory. The consideration of the ZPF 
as a fundamental ingredient of the theory is thus not a subterfuge conceived to 
simplify or guide the calculations, but a fundamental step to unfold the deep 
meaning of the quantum behaviour of matter and field. 

However, it should also be stressed that there can be differences. The present 
theory gives, by construction, only an approximate description of nature. It 
could be that its further development in search of a more detailed or refined 
description leads to discrepancies, open to resolution only by experiment. 
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